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Calculation  of  Threshold  Conditions  for  Materials 
Charging  in  Maxwellian  Plasmas 


1.  BACKGROUND 

DeForest1  found  that  the  ATS-5  spacecraft,  in  geosynchronous  orbit,  charged 
up  to  high  potentials.  Differential  potentials  of  10  to  12  keV  were  found  to  occur 
on  ATS-5  and  ATS-6  in  eclipse,  with  occasional  charging  to  several  kilovolts  in 
sunlight.  DeForest  showed  that  the  theory  of  the  Langmuir  probe  in  a  plasma 
provides  an  explanation  of  the  observed  surface  potentials.  In  contrast  to  labora¬ 
tory  plasma,  geosynchronous  orbit  plasmas  have  temperatures  of  many  kilovolts. 
DeForest  found  that  the  secondary  electron  emission  due  to  these  kilovolt  plasmas 

must  be  included  in  order  to  account  for  the  observed  potentials.  Essentially  the 

2 

same  result  was  found  by  Knott.  Since  that  time,  it  has  been  found  that  the  geo¬ 
synchronous  plasma,  particularly  the  electron  component,  is  typically  dual 
Maxwellian,  having  a  hot  and  a  cold  population,  rather  than  a  single  Maxwellian. 


(Received  for  publication  24  January  1985) 

1.  DeForest,  S.F.  (1972)  Spacecraft  charging  at  synchronous  orbit,  J.  Geophys 

Res.  77;651- 

2.  Knott,  K.  (1972)  Equilibrium  potential  of  a  magnetospheric  satellite  in  an 

eclipse  situation.  Planet.  Space  Sci.  20:1137. 
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Further  charging  of  materials  in  a  dual  Maxwellian  environment  is  a  more  com- 
3  4  5 

plex  phenomenon.  *  '  In  some  cases  two  stable  potentials  are  possible. 

1.1  Charging  in  a  Single- Maxwellian  Flaunt 

In  this  report  we  show  that  materials  in  a  multikilovolt  plasma  have  a 
charging  threshold  in  a  single  Maxwellian  environment.  The  threshold  electron 
temperature  and  charging  curves  are  shown  for  a  variety  of  spacecraft  materials. 
We  show  that  the  controlling  parameter  for  material  charging  in  a  one -Maxwellian 
plasma  is  the  electron  temperature.  In  an  environment  characterized  by  a  single 
temperature,  a  spacecraft  whose  surface  is  composed  of  a  variety  of  materials 
will  charge  differentially.  Materials  with  a  low  threshold  will  be  highly  charged 
while  those  with  a  high  threshold  will  be  uncharged.  In  a  dynamic  situation,  in 
which  the  electron  temperature  is  varying,  because  of  the  existence  of  a  threshold 
and  the  steepness  of  the  charging  curve  above  threshold,  the  potential  of  a  ma¬ 
terial  can  jump  from  zero  to  kilovolts  for  a  small  change  in  electron  temperature. 

1.2  Charging  in  Double-Maxwellian  Plaamaa 

The  dual  Maxwellian  case  is  frequently  encountered  at  geosynchronous  orbit. 
We  derive  the  necessary  and  sufficient  conditions  for  charging  in  a  dual  Maxwell¬ 
ian  environment.  In  order  to  charge  negatively,  the  surface  must  start  out  with 
net  negative  current  on  it  or  it  will  never  reach  a  stable  high  negative  potential. 

In  addition  to  the  electron  temperature  threshold  for  charging  in  a  single  Maxwell¬ 
ian,  there  is  a  threshold  in  the  hot  electron  current  that  must  be  surpassed  in 
order  for  a  surface  to  charge. 


1.3  Charging  Condition* 

In  this  report,  we  consider  a  material  surface  immersed  in  a  plasma  com¬ 
posed  of  two  species;  electrons  and  ions.  The  equilibrium  potential,  V,  of  the 
surface  will  be  determined  by  a  current  balance  equation  of  the  form: 

Je(V)  +  Ji(V)  *  0  (1> 


3.  Prokopenko,  S. ,  and  Laframboise,  J.G.  (1977)  Prediction  of  large  negative 

shade-side  spacecraft  potentials,  Proc.  Spacecraft  Charging  Technology 
Conference,  C.  P.  Pike  andR.R.  Lovell,  Eds.,  AFGL-TR-77-0051  and 
NASA  TMX-73537,  pp.  369-387. 

4.  Sanders,  N.  L. .  and  Inouye,  G.  T.  (1979)  Secondary  emission  effects  on  space 

craft  charging:  energy  distribution  considerations.  Spacecraft  Charging 
Technology  -  1978,  NASA  Conf.  Publ,  2071. 

5.  Besse,  Arthur  L.  (1981)  Unstable  potentials  of  geosynchronous  spacecraft, 

J.  Geophys.  Res.  86:2443. 


where  is  the  current  density  due  to  incident  electrons  and  is  that  due  to  ions. 
These  currents  will  include  contribution  from  secondary  interactions  and  back- 
scatter  in  the  surface  material.  We  will  assume  that  die  initial  voltage  equals 
zero  and  that  the  final  equilibrium  voltage  is  negative.  In  Section  2,  we  write 
down  the  threshold  conditions  for  solutions  to  Eq.  (1)  when  the  plasma  can  be  de¬ 
scribed  by  single  Maxwellian  distributions.  In  Section  3,  we  extend  the  analysis 
to  include  double  Maxwellian  plasmas. 
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Equation  (4)  is  a  Boltzmann  factor  that  reduces  the  flow  of  electrons  because 
of  ihe  repulsive  force.  The  ions  on  the  other  hand  are  attracted  and  the  current 
increases  by  the  orbit-limited  factor  [  Eq.  (5)] . 

The  secondary  and  backscatter  currents  are  proportional  to  the  primary  in¬ 
cident  currents.  Thus  the  net  currents  can  be  expressed  as 

Je(V)  =  jJ(V)  {-1  +  6S(Te)  +  «B(Te)  } 

J^V)  -  JP(V)  {l+fi8^)} 
where  we  define 

g 

6  (T)  ■  electron  secondaries  ratio 
e 

g 

6  (T  )  =  electron  backscatter  ratio 

g 

6  (T.)  =  ion  secondaries  ratio. 

In  Eq.  (7)  we  have  neglected  ion  backscatter,  which  is  a  small  effect.  Also, 

g 

the  ratio  6  (T.)  will  have  a  slight  dependence  on  V,  which  we  can  disregard. 

The  sign  convention  in  Eqs.  (6)  and  (7)  is  that  negative  charge  coming  to  the 
surface  is  counted  as  negative  current.  All  the  6(T)  terms  are  positive  since 
they  correspond  to  cases  where  electrons  are  ejected  from  the  material. 

Substituting  Eq.  (2)  into  Eq.  (6)  yields: 

/  kT 

Je W  m  Ne  (a TZTJ  ql'1  +  6  (Te>  +  ‘X^e^  <8) 

and  putting  Eq.  (3)  into  Eq.  (7)  gives 
I  kT 

Jl(V>  -  Nl [ukj  11  •  g,<V>  .  (9) 

We  shall  assume,  in  accordance  with  Eqs.  (4)  and  (5),  that  g  (V)  and  g^V) 
are  positive  functions  and  that  gg(V)  approaches  zero  monotonically  as  -V  —  », 
but  gj(V)  is  non-decreasing.  Our  results  will  thus  be  valid  for  any  probe  geometry 
for  which  this  is  true  (such  as,  spherical,  cylindrical,  or  planar). 

It  is  evident  from  Eq.  (9)  that  Jj(V)  is  always  positive  since  this  is  true  for 
fj(V)  and  for  the  bracket  factor.  In  order  for  overall  current  balance  to  occur  in 
Eq.  (1),  we  must  have  Je(V)  <  0.  But  fi'om  Eq.  (8),  the  only  way  for  this  to  occur 
is  if  the  condition 


(6) 

(7) 


i 


«(Te)  =  -1  +  6S(Te)  +  6B(Te)  (10) 

is  satisfied.  The  function  6(Te)  is  plotted  in  Figure  1  for  five  different  materials. 
For  all  temperatures  below  the  point  T*  where  6(Tg)  crosses  zero,  no  charging 
is  possible.  The  symbols  and  material  properties  assumed  are  listed  in  Table  1. 


T«  (K«V) 

Figure  1.  Charging  Thresholds  for  Materials  in  a  One- 
Maxwellian  Plasma.  The  charging  threshold  is  that 
electron  temperature  at  which  the  threshold  factor  is 
zero 

The  terms  6®(Te>  and  6B(Te)  are  shown  separately  in  Figures  2  and  3.  The  sec¬ 
ondary  interaction  terms  6s(Te)  become  large  at  low  temperatures  while  the  back- 
scatter  term  6  (T g )  is  relatively  flat.  When  Tg  <  Tg*  sum  Qf  these  coeffi¬ 
cients  becomes  greater  than  1  and  thus  more  electrons  are  ejected  from  the 
surface  than  are  arriving  from  the  plasma.  The  net  effect  in  this  case  is  to 
discharge  the  surface. 

Equation  (10)  is  a  necessary  but  not  a  sufficient  condition  for  charging.  Jg(V) 
must  not  only  be  negative  but  it  also  must  be  large  enough  to  cancel  J.(V).  From 
our  assumption  that  gg(V)  —  0  as  -V  -*  »  it  follows  from  Eq.  (8)  that  |<Je(V)| 
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Table  1.  Material  Parameters 


Material 

Plot 

Symbol 

Electron 

Secondary 

6  E 
m  m 

Ion 

Secondary 

1  keV  proton 
yield  (dE/dX)m 

Z 

Kapton 

o 

2.1.  0.15 

0.455,  140 

5 

SOLAR 

A 

2.05,  0.41 

0.244,  230 

10 

Astroquartz 

+ 

2.40,  0.40 

0.455,  140 

10 

Aluminum 

X 

0.97,  0.30 

0.244,  230 

13 

Gold 

0 

0.  88,  0.  80 

0.43,  135 

79 

SOLAR  refers  to  solar  cell  cover  glass  material 

Astroquartz  is  a  quartz  cloth  material 

1.0  2.0  3.0  4.0  S.O  6.0  7.0  8.0  9.0  10.0 


Te  {KeV) 

Figure  2.  Electron  Secondary  Emission  Yields  for  Five 
Materials 


0.900 


o  KAPTON 


Figure  3.  Electron  Backscattering  Coefficient  for  Five 
Materials.  Astroquartz  and  SOLAR  are  coincident 


On  the  other  hand,  we  have  seen  that  ,L(V)  is  always  positive.  Therefore  a  neces¬ 
sary  and  sufficient  condition  for  charging  is: 

-Je<0)  >  J.(0)  (11) 

for  then,  as  -V  increases,  there  will  be  some  point  at  which  the  function  -J  (V) 
crosses  over  J.(V),  that  is,  the  current  balance  in  Eq.  (1)  will  be  satisfied. 

From  Eqs.  (8)  and  (9)  we  can  write  the  conditions  of  Eq.  (11)  as 

m 

-N  T  6(T  )  >  ~~  N.T.6(T.)  ,  (12) 

e  e  e  rm  i  t  t 

where  «(T.)  =  1  +  6S(T.). 

The  ion  secondaries  coefficient  6®(T^)  is  given  in  Figure  4  for  the  five  mater¬ 
ials  The  effect  of  6(T.)  is  reduced  because  of  the  multiplicative  factor  mg/mi 
((183R)'^2  for  protons] .  If  the  ion  term  is  neglected  we  recover  the  previous 
condition  (  Eq  ( 10)]  . 

The  dynamics  of  the  charging  is  illustrated  in  Figure  5,  which  shows  -J  (V) 
and  L(V)  versus  V  for  astroquartz  material.  Each  curve  corresponds  to  a  differ¬ 
ent  Maxwellian  temperature  (T  ,  T.  =  5,  10,  15,  20,  2  5).  The  plasma  densities 


ELECTRON  AND  ION  CURRENTS 


Figure  4.  Ion  Secondary  Emission  Yields  for  Five  Materials 


-2.5  -  5.0  -7.5  -10.0-12. 5-15.0-17. 5-20  0-22. 5-25.0 
V  (  K  V  ) 


Figure  5.  Electron  and  Ion  Currents,  Single  Maxwellian 
The  equilibrium  surface  voltage  is  defined  by  the  inter¬ 
section  of  the  electron  and  ion  current  curves 


i 


i 


are  taken  as  Ne  =  N^  =  1.  The  intersection  of  any  two  curves  determines  the 
equilibrium  voltage  V.  The  curve  for  -Je(V)  with  Tg  «  5  lies  below  the  origin  and 
hence  can  never  intercept  a  J^V)  line.  This  is  an  example  of  a  case  where  the 
electron  temperature  is  below  the  threshold,  that  is,  Eq.  (10)  is  not  satisfied  and 
no  charging  can  occur.  The  remaining  -Je(V)  curves  start  above  the  Jj(V)  lines  at 
V  ■  0  so  that  the  condition  of  Eq.  (11)  holds  and  charging  takes  place. 

Note  that  there  is  a  greater  separation  (AJ/&T)  of  the  electron  curves  than 
there  is  for  the  ions  and  therefore  the  electron  temperatures  have  a  stronger  in¬ 
fluence  on  the  equilibrium  voltages .  This  comes  about  because  the  ions  have  much 
larger  mass  and  thus  the  ion  currents  are  smaller  by  a  factor  me/rm  relative  to 
the  electron  currents  [see  Eqs.  (8)  and  (9)].  These  computations  were  carried 

out  with  a  materials  charging  code,  MATCHG,  which  is  a  subset  of  the  NASCAP 
0 

vehicle  charging  code. 

3.  THRESHOLD  CONDITIONS  FOR  DOUBLE  MAXWELLIAN  PLASMAS 

For  a  multi -Maxwellian  plasma,  Eqs.  (8)  and  (9)  become: 
n 

Je(V)  =  £  N3e  (kT3e/2Tme)l/2  6(T3e)gJe(V) 

J-l 

and 

Q 

Ji(V)  =  Z  Nj(kT{/2ffm.)1/'2  6(T|)g{(V) 

J-l 

where  the  superscript  j  labels  the  contribution  from  each  Maxwellian  component. 
As  before,  the  ion  current  J.(V)  is  always  positive  since  each  term  in  the  summa¬ 
tion  is  greater  than  zero.  Again  we  must  have  Jg(V)  <  0  for  charging  to  occur. 

The  necessary  condition  for  charging  is  more  complicated  than  for  a  single  Max¬ 
wellian.  At  least  one  of  the  T3  must  be  above  threshold  and  the  sum  of  all  the 

e 

contributions  to  Jg(V)  must  be  negative. 

We  can  see  that  J  (V)  -»  0  as  -V  —  «  because  each  ge(V)  —  0  and  the  threshold 
conditions  for  multi-Maxwellian  can  be  again  taken  to  be  Eq.  (11).  Writing  out 
Eq.  (11),  using  Eqs.  (13)  and  (14),  we  get 

6.  Katz,  I. ,  Cassidy,  J.  J. ,  Mandell,  M.  J. ,  Schnuelle,  G.  W. ,  Steen,  P.  G. , 
and  Roche,  J.  C.  (1979)  The  capabilities  of  the  NASA  charging  analyzer 
program.  Spacecraft  Charging  Technology  -  1978,  NASA  Conf.  Publ.  2071, 
AFGL-TR  -75082! 
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(13) 


(14) 


(15) 


p  Q 

-  Z  nJ  (tJ)1^2  6  (tJ)  >  (me/m.)l/2  Z  n{(t{)1/2  6  (t{) 
i=i  j=i 


For  P  =  Q  =  1  we  recover  the  previous  result  [Eq.  <  12)] .  If  the  electron  distri¬ 
bution  can  be  described  by  a  double  Maxwellian  (m  =  2)  then  we  find 

N JTy/2  «(Tg)  +  N2(T2)1//2  6  (T2)  <  (m  /m  J1/2  F.  .  (16) 

C  C  C  C  C  C  v  1  1 

where  we  have  defined 


F.  = 


£  n|  6  (T^)  (Tj)1^2 


(17) 


j=l 


The  quantity  F.  is  always  negative.  Thus  to  satisfy  Eq.  (16),  at  least  one  of  the 
1  (21 

electron  temperatures,  say  T  must  be  above  the  threshold.  From  Eq.  (16), 


the  borderline  in  (N 


(2)  M(l) 


N' *')  space  that  separates  charging  and  non-charging 


cases  can  be  written  as: 


+  B 


with  slope  and  intercept  given  by: 


(T<“) 

Il/2  1 
1  6  ( 

>«») 
t  e  / 

1/2 

6 

U2)\ 

V  e  ) 

(  me  > 

V  IMj  ) 

f 

1/2 

1  Fi 

\  e  J 

(18) 


(19) 


(20) 


4.  DISCUSSION 

In  the  dual  Maxwellian  situation,  a  typical  case  is  selected  to  show  threshold 
effects.  The  material  is  SOLAR,  or  solar  cell  material.  A  single  ion  tempera¬ 
ture  of  10  keV  is  employed.  A  low  temperature  electron  Maxwellian  at  1  keV  is 
chosen.  We  show  the  results  of  variations  of  temperature  and  densities  of  the  hot 
electron  component,  since  that  is  the  critical  component  for  charging. 
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Figure  6  shows  the  boundary  lines  in  (N^2A  N^)  space,  defined  by  Eq.  (18), 

/i\  i  e  (2) 

for  aluminum.  T  has  been  chosen  below  the  threshold  temperature  and  T  is 
above  it.  The  regions  above  the  straight  line  contain  those  density  points  for 
which  charging  occurs. 


Figure  6.  The  Boundaries  in  (p/2\  N^)  Space  Separating  Charging  and 
Non -Charging  Conditions  in  a  Two -Maxwellian  Plasma 

Figure  7  shows  the  current  voltage  curves  for  a  single  ion  temperature  and 

four  cases  of  two -Maxwellian  electrons.  The  low  electron  temperature  has  been 

chosen  to  be  1  keV  below  the  charging  threshold.  The  high  electron  temperatures, 
(2) 

T  ,  are  5,  10,  15  or  20  keV.  The  charging  potentials  correspond  to  the  points 
at  which  the  electron  and  ion  currents  intersect,  since  this  is  a  graphical  version 
of  the  condition  that  the  net  current  is  zero. 

There  are  four  cases.  For  Tg  =  1,  20  there  is  a  single  root  at  negative  po¬ 
tential.  The  T  =  1,  15  case  has  two  roots,  one  stable  at  V„  and  an  unstable  root 
at  Vg-  In  the  Tfi  =  1,  10  case  the  2  roots  at  Vg  and  Vg  are  close  together,  and 
nearly  coalesce.  The  fourth  case  is  for  T  ■  1,  5  where  there  are  no  roots  at  all. 

Starting  at  zero  potential  on  the  material,  only  one  case  will  be  negatively 
charging.  This  is  the  case  in  which  for  zero  surface  potential  the  electron  cur¬ 
rent  exceeds  the  ion  current;  otherwise,  the  surface  would  charge  positively.  In 
a  two -Maxwellian  plasma,  not  only  must  there  be  a  current  balance  for  negative 
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V  (Kilovolt*) 


Figure  7.  SOLAR  Material  in  a  Two -Maxwellian  Plasma. 

The  equilibrium  surface  voltage  is  again  defined  by  the  in¬ 
tersection  of  the  electron  and  ion  current  curves.  See  the 
text  for  further  details 

charging  to  occur,  but  also  the  surface  current  must  start  out  with  net  negative 
current  or  it  will  never  reach  the  stable  high  negative  potential. 

5.  CONCLUSIONS 

We  have  examined  the  threshold  conditions  for  charging  of  a  surface  material 
immersed  in  a  Maxwellian  plasma.  The  threshold  relations  depend  on  the  plasma 

temperatures  and  densities  and  on  the  material  properties  through  the  secondary 

S  S  B 

emission  coefficients  6  (Te),  6  (T.)  and  the  backscatter  ratio  6  (Te).  For  a 

single  Maxwellian  plasma,  the  necessary  condition  for  charging  is  that  the  elec¬ 
tron  temperature  must  be  greater  than  the  threshold  temperature  defined  by 

6  (Tg)  =  -1  +  6S(Tg)  +  6B(Te)  =  0 

For  a  multi-Maxwellian  plasma  at  least  one  of  the  electron  temperatures  must  be 
above  threshold  and  the  net  electron  induced  current  must  be  negative.  To  arrive 
at  necessary  and  sufficient  conditions  for  charging,  one  has  to  take  into  account 
the  ions  as  well  as  the  electrons.  Explicit  formulas  have  been  given  for  single 
and  double  Maxwellians  [see  Eqs.  (12)  and  (16)].  For  a  double  Maxwellian,  the 
threshold  condition  defines  a  linearly  bounded  region  in  the  density  space 
(N(£>,  N(£}). 
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